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Abstract 

This  paper  presents  an  efficient  algorithm  for  solving  a class  of 
quadratic  Integer  programming  problems.  These  problems  Include  discrete 
versions  of  the  quadratic  placement  problem  and  the  squared  Euclidean 
distance  problem.  The  algorithm  solves  a finite  sequence  of  minimum  cut 
problems,  or  equivalently  maximum  flow  problems,  on  a graph  with  n + 2 
vertices  where  n Is  the  number  of  variables  In  the  problem. 


Consider  an  Integer  program  of  the  form 


U 11  11 

minimize  f(Y)  - I q . y y.  + I b y 

j-1  k-1  J * J-1  ^ J 


s.t.  

where  Aj  and  are  non-negative  integers,  - ‘^kj  ’ **jk  — ^ j ^ 
n 

J q , ^ 0 for  j ■■  1,  2 n.  We  will  develop  an  algorithm  which  solves  (PO) 

J-1  ^ 

n 

by  solving  a sequence  of  not  more  than  J (u.  “ minimum  cut  problems  on  a 

j*l  ^ ^ 

graph  with  n+2  vertices. 

Problems  which  are  of  Che  form  (PO)  include 
discrete  versions  of  the  squared  Euclidean  distance  location  problem  (White 
[6])  and  the  n-dimenslonal  quadratic  placement  problem  (Hall  [1]).  The  dis- 
crete squared  Euclidean  distance  location  problem  has  the  form 

j*l  k»l  J J J j«l  k=l  -J  J J 

..t.  Zj  . tj.  tjtl.  ....  Uj  znd  zj  - u’ 

I * 

where  (zj,  Zj)  is  the  location  of  new  facility  j,  (a^,  a^)  is  the  location  of 
old  facility  1,  and  v..  and  w are  non-negative  weights  normally  corresponding 
to  Che  number  of  crips  between  facilities.  This  problem  decomposes  into  a 

t 

problem  in  z and  a problem  in  z both  of  which  are  easily  put  into  the  form  (PO) . 
The  n-dlmensional  quadratic  placement  problem  can  be  considered  as  a variant 
of  the  squared  Euclidean  distance  problem. 

The  underlying  ideas  for  Che  algorithm  Co  be  presented  here  are  as 


1 


follows.  We  note  from  Picard'  and  Ratliff  [2]  that  a restricted  version  of 


(PO) , where  each  can  take  on  only  its  lower  bound  or  its  lower 

bound  plus  one,  can  be  efficiently  solved  as  a minimum  cut  problem  on  a graph. 
(This  will  be  Illustrated  later.)  If  Y*  is  the  optimum  solution  to  this  re- 
stricted problem,  we  will  show  that  there  exists  an  optimum  solution  Y®  to 
the  original  problem  with  y°  ^ Yj  2,  ...,  n.  Hence  for  each  j with 

y*  » Aj+1  we  can  Increase  the  lower  bound  on  y^  by  one.  We  can  then  repeat 
this  process  until  we  obtain  an  optimum  solution  Y*  with  each  y*  equal  to  its 
lower  bound.  We  will  show  that  such  a solution  is  optimum  to  the  original 
problem.  Before  formalizing  the  algorithm,  we  will  develop  the  necessary 
theoretical  underpinning. 


Theoretical  Results 

Consider  the  problem 
minimize  f(X) 

s.t.  Xj  = 0,  1,  2 uj  for  j * 1,  2,  ...,  n 

where  f(*)  is  as  defined  for  (PO)  and  uj  is  a nonnegative  Integer  for 
J » 1,  2,  ...,  n.  Also  consider  the  related  problem 
minimize  f(X) 

s.t.  Xj  = 0 or  1 for  j = 1,  2,  ...,  n 
where  f(«)  is  as  defined  for  (PO) . 


(PI) 


(P2) 


Lenrna  1;  If  X*  is  an  optimum  solution  to  (P2) , then  there  is  an  optimum  solu- 
tion X'  to  (PI)  with  Xj  i 

Proof  of  Lemma  1:  Let  X®  be  any  optimum  solution  to  (PI)  and  let  R ■ t J l*j  * 1 
and  xj  - 0}.  Define  the  vectors  X*  and  X"  as 
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1 if  j £ R 

X®  otherwise. 
‘ j 

0 if  j e R 


X*  otherwise. 


Noting  that  for  each  jfR  the  variables  x*  = 1 and  xj  » 0 we  have 

f(X*)  - f (X")  ~ I I q +2  ^ I "Ilk  I b < 0 

jeR  k£R  j£R  k^R  ^ j£R  ^ 


(Rl) 


since  X"  is  feasible  to  (P2)  and  X*  is  optimum  to  (P2) . Also  noting  that  for 
each  jeR  the  variables  x®  = 0 and  Xj  = 1 we  have 

(R2) 


f(X’)  - f(X®)  = I Iq.^+2l  I q^k  *k  ^ '’l 
j£R  k€R  ^ jeR  k^R  ^ j€R  ^ 


Now 


since  q , £ 0 and  for  k^  we  have  x®  ^ x*  it  follows  that  each  term  of 

J J 

(R2)  is  less  than  or  equal  to  the  corresponding  term  of  (Rl) . Hence, 

f(X')  _<  f(X®).  Since  X'  is  feasible  to  (PI)  and  has  xj  ^ x*  the  result  follows. 

Q.E.D. 

Lemma  2;  If  X*  is  an  optimum  solution  to  (P2)  and  x*  * 0 for  j=l,  2,  ...,  n, 
then  X*  is  optimum  to  (PI) . 

Proof  of  Lemma  2:  Assume  that  there  exists  an  X®  feasible  to  (PI)  with 
f(X®)  < 0.  (Note  that  f(X*)  - 0.) 

Define  the  vectors  X'  and  X"  as 


x'  = < 

j ' 


1 if  xj  > 1 
0 if  xj  * 0 


xV  » X?  - x!  for  j ••  1,  2,  ...» 


j j j 


Then 


n n 


f(x®)  - f(x’  + X")  ~ I I q..  (x'+x")(x;+x;;)  + I b (x'+x") 

j-1  k-1  J J J j-1  ^ J J 

n n 

C)  + f(X")  + 2 ^ J q x’  xV  < 0 

j-1  k-1  J 


f(X' 


(R3) 
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Since  X'  is  feasible  to  (P2)  and  X*  is  optimum  to  (P2)  with  f(X*)  ■ 0,  it 
must  be  the  case  that  f(X')  > 0. 


n n 

Now  consider  the  term  2 ^ ^ q , x'  x". 

J-1  k-1  J J 


From  the  assumption  that 


n n 

1 ^ 0 we  have  ^ q . x'  ^ 0 if  xl  ■ 1.  In  addition,  whenever  x*  ■ 0 

k-1  k-1  J J J 


n n 

we  have  x"  “ 0.  Hence,  2 ^ ^ q , x'  xV  ^ 0.  Therefore,  from  (R3)  it  follows 

J J-1  k-1  J 

that  f(X")  < 0. 

If  we  repeat  this  process,  each  time  setting  X**  - X",  after  a finite  num- 
ber of  repetitions  we  must  obtain  an  X"  with  f(X")  < 0 and  x'J  = 0 or  1.  This 
contradicts  the  assumption  that  X*  is  optimum  to  (P2) . 

Q.E.D. 

We  need  to  make  one  further  observation  before  stating  an  algorithm. 

If  a problem  is  in  the  form  (PO) , we  can  make  the  change  of  variables 
Xj  - to  obtain 

n n n n 

min  f(X)  + L - ^ [ q x x.  + J]  x (b  + 2 J] 

j-1  k-1  J j.l  J J k-1  J 

s.t.  Xj  - 0,  1,  ...,  Uj-fj. 

Since  the  constant  term  does  not  affect  the  optimization,  (P3)  is  of  the  same 
form  as  (PI) . Note  that  the  only  coefficients  of  (P3)  which  differ  from  (PO) 
are  those  for  the  linear  terms. 


Algorithm 

(1)  Given  a problem  of  the  form  (PO) , transform  it  via  the  change  of  variables 
Xj  - y^  - f j to  the  form  (PI) . 


A 


(2)  Solve  the  corresponding  problem  (P2)  as  a minimum  cut  problem  on  a graph 
(e.g.,  Picard  and  Ratliff  [4])  to  obtain  the  solution  X*.  Let 


S - {jIx"  - 1}. 


(3)  If  X 


stop.  It  follows  from  Lemma  2 and  the  transformation  used  to 


obtain  (P3)  that  for  j ■ 1,  2,  . . . , n Is  optimum  to  (PO) . 


(4)  If  S set  Jlj  ■ 2.J  + 1 for  all  j e S.  It  follows  from  Lemma  1 and 


the  transformation  used  to  obtain  (P3)  that  the  new  is  a valid  lower 


bound  for  In  (PO) . For  each  j such  that  Uj  permanently  set 


yj  “ . Go  to  step  (1) , 


Since  at  least  one  lower  bound  la  Increased  by  one  at  each  iteration,  the 


algorithm  terminates  after  at  most  ^ (u  -1  ) iterations. 

j=l  ^ ^ 


Example  Problem 

In  order  to  demonstrate  the  steps  of  the  algorithm,  consider  the  following 
example  problem: 


~6  -4  -T 

■-6" 

minimize  [y^^,  y2,  y3] 

-4  8 -4 

^2 

+ [yj^.  72*  73! 

-8 

-2  -4  20 

-4 

— 

J ^ 

• * 

s.t.  - 0,  1,  2,  3 y2  “ 0, 

1. 

2 

y^  - 0,  1,  2,  3, 

4 

(P4) 


Note  that  this  problem  Is  already  In  the  form  (PI) , so  no  Initial  change  of 
variables  Is  required. 

In  order  to  solve  the  corresponding  problem  (P2)  we  can  define  an  undirected 
graph  with  vertices,  0,  1,  2,  ...,  n+1  and  capacities  defined  as  follows 
(Picard  and  Ratliff  [5]). 
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• -‘‘IJ 


il 

=J.n+l  ■ "jk  ^ "j-  “> 


'oj  * “’‘‘-Ji''dk  - "j*  '» 

If  (S,  S)  Is  a minimum  capacity  cut  in  the  graph  separating  node  0 and  node  n-fl 
with  0 £ Sand  n+1  € S,  an  optimum  solution  to  (P2)  is  • 1 for  j £ S and 
Xi  - 0 for  j £ S. 

For  this  example  the  graph  Is  as  shown  In  Figure  1 where  the  dashed  lined 

Indicates  a minimum  cut.  The  optimum  solution  to  (P2)  Is  x^  > 1,  x^  * 1,  and 

x^  * 0.  From  step  4 of  the  algorithm  we  set  ■ 1,  ^ ■ 0. 

Performing  the  change  of  variables  (P3)  on  (PO)  with  the  new  lower  bounds 

T 

yields  the  new  linear  term  coefficients  [-2,  0,  -16]  . The  constant  term  Is 
Ignored  and  the  remaining  coefficients  are  the  same  as  In  (P4).  The  new  graph 

and  minimum  cut  are  shown  In  Figure  2.  An  optimum  solution  to  the  new  problem 

(P2)  Is  Xj^  “ 1,  X2  ■ 1,  and  x^  “ 1.  From  step  4 we  get  * 2,  2,  and 

= 1. 

Again  performing  the  change  of  variables  (P3)  on  (PO)  with  the  new  lower 
bounds  yields  the  new  linear  term  coefficients  [-2,  0,  12]  . The  new  graph 
and  minimum  cut  are  shown  In  Figure  3.  This  time  the  optimum  solution  to  (P2) 

Is  Xj^  - 0,  X2  * 0,  and  x^  ■ 0.  Hence,  an  optimum  solution  to  (PO)  is  “ 2, 

y2  “ ^2  “ ^3  " ^3  “ 


1 


min  cut 

/ 


Figure  2:  Graph  for  example  Iteration  two. 


Figure  3:  Graph  for  example  Iteration  three. 


Constrained  Problems 

Constraints  of  the  form  yj  ~ Integer  can  be  Incorporated 

Into  (PO)  very  simply  by  using  a device  developed  In  Hammer,  Rosenberg,  and 

Rudeanu  [2],  and  discussed  In  Hammer  and  Rudeanu  [3,  p.  123].  For  each  con- 

+ **  2 

stralnt  of  this  form  we  add  to  f(y)  In  (PO)  the  term  (S  - S 

where  and  S are  the  sums  of  the  positive  and  negative  coefficients  of  f(y), 

respectively.  If  the  optimum  solution  to  the  new  problem  Is  feasible  to  the 

constraints.  It  Is  optimum  to  the  original  problem.  If  the  optimum  solution 

to  the  new  problem  Is  not  feasible  to  the  constraints,  then  there  does  not  exist 

a feasible  solution  to  the  original  problem.  Note  that  (S^  - S + 1)  Is  posl- 

+ - 2 

tlve,  hence  adding  terms  of  the  form  (S  - S + 1)  (Yj  “ Yfc  " 
yields  a function  with  the  same  form  as  the  original  f(y). 

By  Including  constraints  of  this  form  In  the  squared  Euclidean  location 
problem,  we  are  able  to  specify  the  distance  between  pairs  of  facilities.  This 
device  can  be  used  to  model  the  location  of  line  segments  (rectangles  for  the 
two  dimensional  problem)  by  considering  each  end  of  the  segment  as  a point  to 
be  located  and  then  constraining  the  distance  between  these  points  to  be  the 
length  of  the  segment.  This  Is  of  particular  Interest  In  modeling  layout 
problems. 


Conclusions  and  Extensions 

The  algorithm  presented  here  can  be  considered  as  a generallzatlor  of  the 
algorithm  In  Picard  and  Ratliff  [5]  for  the  rectilinear  distance  location  pro- 
blem. The  basic  difference  Is  that  additional  special  structure  In  the  recti- 
linear problem  (l.e.,  each  cut  problem  solved  partitions  the  variables  Into 
two  sets  and  each  set  can  be  solved  as  an  Independent  problem)  allows  a more 
efficient  algorithm  which  requires  at  most,  m-1,  where  m Is  the  number  of  old 
facilities,  minimum  cut  problems  on  a graph  with  n+2  vertices. 
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